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ABSTRACT 

We study the baryon content of low mass galaxies selected from the Sloan Digital Sky Survey 
(SDSS DR8), focusing on galaxies in isolated environments where the complicating physics of galaxy- 
galaxy interactions are minimized. We measure neutral hydrogen (HI) gas masses and line-widths 
for 148 isolated galaxies with stellar mass between 10^ and lO^'^M©. We compare isolated low mass 
galaxies to more massive galaxies and galaxies in denser environments by remeasuring HI emission 
lines from the Arecibo Legacy Fast ALFA (ALFALFA) survey 40% data release. All isolated low mass 
galaxies either have large atomic gas fractions or large atomic gas fractions cannot be ruled out via 
their upper limits. We measure a median atomic gas fraction of /gas = 0.81 ± 0.13 for our isolated 
low mass sample with no systems below 0.30. At all stellar masses, the correlations between galaxy 
radius, baryonic mass and velocity width are not significantly affected by environment. Finally, we 
estimate a median baryon to total dynamical mass fraction of /baryon,disk = 0.15 ± 0.17. We also 
estimate two different median baryon to halo mass fractions using the results of semi-analytic models 
(/baryon,halo = 0.04 ± 0.06) and abundance matching (/baryon,halo = 0.04 ± 0.02). Baryon fractions 
estimated directly using HI observations appear independent of environment and maximum circular 
velocity, while baryon fractions estimated using abundance matching show a significant depletion of 
baryons at low maximum circular velocities. 

Subject headings: galaxies: dwarf - galaxies: evolution - galaxies: ISM - galaxies: kinematics and 
dynamics 


1. INTRODUCTION 

Environment plays a major in role in th e evolution and 
the observed properties of galaxies (e.g.. Dressier 1980 


baryons in atomic gas (e.g., Grcevich & Putman 2009 


et al. 


Kauffmann et an2004[ Blanton et al. ___ 

2010 1 |Geha et ar.||2012[ |Tal et al.||2bl4[ IvVetzel et al. 

2014k . Environmental processes such as tidal forces and 


ram pres sure stripping shape , quench and destroy galax- 


ies (e.g., Toomre & Toomre 

1972 Gunn & Gott 

1972 

Moore et al. 1996 Beharrubia et al. 

2UU8). Thes 

s pro- 


with stellar masses below 10^ Mq, where more massive 
nearby galaxies dramatically affect the local gravitational 
potential (e.g. , Kra vtsov et al.||2004 Zolotov et al.||2012 
Kenney et al.||2013l). 


'The effects of environment complicate both observa¬ 
tional interpretations and comparisons to theoretical pre¬ 
dictions of galaxy formation. Galaxies in isolation are rel¬ 
atively undisturbed by environmental processes. Thus, 
isolated galaxies can pote ntially present a clearer pic- 


B art on 


et al. 120071 ICortese et al. 120111 IKarachentsev et al.|2Ullt 

Tollerud et al. 2011 Toribio et al. 

2011|). Recent work 

has found that the properties of iso 
ies, loosely defined as “field galaxie 
with small-scale predictions of ACE 

lafed low mass galax- 
s”, appear to conflict 
)M cosmology (IKirby 

et al. 

2014 

Klypin et al. 

2014 

Papastergis et al. 

2015 . 


defined and well studied. 

Isol ated galaxies retain t heir HI gas out to larger radii 
(e.g.. Begum et al. 20081, contain a large fraction of 


Geha et al. 2006 hereafter G06) and are predicted to 


exhibit larger maximum circular velo cities than galaxies 
that have been tidally stripped (e.g. 

2014|. Therefore in underdense regions 


Brooks & Zolotov 
where galax¬ 
ies are less affected by environmental processes, it might 
be easier to predict their properti es when compared 
to galaxies in dense environments (McGaugh fc Wolq 
20101 [Cortese et al.||2 0111 IHuang et al.||2012a| lArgudo- 
Ferndndez et al.||2014p7 


Here we present the largest homogeneously measured 
data set of baryons in isolated low mass galaxies with 
IO'^Mq < M* < 1O® ®M 0. Our new data probe isolated 
galaxies with average circular velocities of 45 km s“^ and 
average baryon masses of 2 x 10® Mq, with circular ve¬ 
locities as low as ~ 20 km s“^ and baryon masses as low 
as ~ IO^'^Mq. We complement th e existing ALFALEA 
40% data release of Haynes et al. (2011), hereafter AL- 
FALFA40, with new 21 cm observations, incorporating 
our isolation criterion into the selection of 148 low mass 
targets. 

This paper is organized as follows: In §[^ we introduce 
the data, define our isolation criteria and describe our au¬ 
tomated HI parameter estimates. In §[^we present our 
results, examining the relationship between galax y ste llar 
mas s, atomic gas mass and atomic gas fraction (§ 3.1 and 


3.21, the relationship b etwe en galaxy size, baryonic mass 
and dyn amical mass (§3.3), and finally baryon fractions 
in § |3.4| In §i we discuss our results in context with the 
literature, particularly with respect to the effects of en- 
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vironment on galaxy properties. Finally, we summarize 
our findings in §[^ 

Throughout this work we assume the following cosmo¬ 
logical parameters: flo = 0.3, = 0.7, Hq = 70 km s“^ 

(i.e., h = 0.7), and a co smic baryon fraction of ft — 0.17 
(Komatsu et al. 2011). We present a small sample of 
isolated low mass galaxies for reference in Table [J Since 
our HI observations are ongoing, our complete data set 
and raw HI spectra are available upon request. This pa¬ 
per is the first in a series using isolated low mass galaxies 
as a control sample for testing galaxy evolution at small 
scales. 

2. DATA 

2.1. Galaxy Catalog 

To construct a sample of isolated low mass galaxies, 
we begin with the NASA-Sloan Atla^ (NSA) catalog 
(|Blanton et al.|2011|), which is itself a re-reduction of the 
SUbb Data Kelease 8 (^DR8, Aihara et al.|20 lll). This 
sample was first described in (deha et al. 15072 Tiereafter 
G12), and we provide a brief description here. 

The NSA re-reduction is opti mized for nearby ex - 
tended objects with z < 0.055 (Blanton et al. 2011). 
The NSA photometry includes image mosaics from both 


is common to define an environment var iable, such as 
nth n earest neighbor surface density (e.g., Brough et al . 
2013) or tidal force index (e.g., Karachentsev et al.|201ip . 


SPSS ugriz and GALEX FUV, NUV bands (Martin et al. 
2005). Sersic model structural parameters are ht using 


the SPSS r-band image. The SPSS spectra are also rean - 
alyzed and recalibrated ( Yan|2011 Yan & Blanton 2012). 
The spectra equivalent widths are re-measured and in- 
cluded in the NSA catalog. 

In order to ensure accurate distance estimates and to 
stay within the spectroscopic completeness limits of the 
SPSS, we discard any galaxy in the NSA catalog where 
it is not detected in the SPSS g, r or i passbands and 
only include galaxies with spectroscopy from the main 
SPSS legacy survey with z > 0.002 and < 17.72. We 
also discard galaxies where the Sersic half-light radius 
is larger than the Petrosian 90% light radius, which in¬ 
dicates that surface brightness model fitting has failed. 
Our final catalog is comprised of 106374 galaxies that 
span distances of 10 < P < 250 Mpc and stellar masses 
of IO^ ^Mq < M* < 1O“-9M0. 

Galaxy stellar masses (M*) are calc ulated in the NSA 
catalo g using the kcorrect software of |Blanton fc Roweis| 
(|2007|) using t he SPSS a.nd GA LEX photometric bands 
and assuming a Ghabrier (2003) IMF. The error in stellar 
mass is estimated from the scatter between models ex¬ 


amined by Blanton & Roweis (2007) who found roughly 
uniform differences of 0.2 dex between models (see their 
Fig. 17). We calculate luminosity distances using SPSS 
re dshifts corrected for peculiar velocities using the model 
of|Willick et al.| (|I997|) , who compiled a velocity-distance 
catalog of spiral galaxies using Tully-Fisher calibrated 
distances and redshift velocities. Pistance uncertainties 
are folded into all distant dependent calculations. 

2.2. Galaxy Environments, Isolation Criteria and 
Sample Definition 

Many methods have been devised for quantifying the 
local and large-scale environments around galaxies. To 
measure environment density and degree of isolation, it 

^ http://www.nsatlas.org 


It is also common to define galaxies “in the held” in or- 
der to study galaxies in isolation, where such a galaxy 
has various definitions in the literature. Other studies 
measure the density of galaxies within a fixed-aperture 
on the sky or measure 2P projected distances bet ween 
galaxies with cuts applied in redshift space (e.g., Gro- 
|ton et al.||2005[ IConroy et alj 2(3()7| [Gallazzi et M. |%U8[ 

3e IMuldrew et 


G12). See Muldrew et al. (|201I|) for a comparison of 

20 different environment definitions applied to a mock 
galaxy catalogue. 

Our goal is not to review and test the various methods 
of quantifying environment. Our goal is to select the 
most isolated galaxies where environmental processes are 
minimized. GI2 has found that low mass galaxies with 
IO^Mq < M* < 10^Mq are always star-forming when 
located more than 1.5 Mpc from a massive host. In our 
work, we assume that environment has not significantly 
affected such low mass galaxies and we adopt the results 
of G12 as our isolation criter ion for low mass galaxies 
(also see Wheeler et al. (2014)). 


For low mass galaxies with < 10^'^Mq, we cal¬ 
culate the 2P projected distance (dhost) to the nearest 
more massive host within a redshift of 1000 km s“^ and 
within a 2P projected distance of 7 Mpc. Massive neigh¬ 
bors may be located beyond the edges of the SPSS sur¬ 
vey, therefore we also search for massive hosts beyond 
the SPS S footprint using the 2MASS Extended Source 
Gatalog ( [Jarrett et al. 2000). A galaxy is considered 
massive if M* > 2.5 x IO^'^Mq (the 2MASS magnitude 
< —23). Therefore, massive hosts are at least ~ 0.5 
dex more massive than the low mass galaxies. Redshifts 
for 2MASS hosts are obtained using a variety of other 
surveys and literature as described in G12. We therefore 
define a low mass galaxy to be isolated if dhost >1-5 Mpc. 
We set dhost = 7 Mpc if no such nearby massive host 
is found. 

Since the results of G12 apply only to low mass galax¬ 
ies, the above isolation criterion is ill-defined for massive 
galaxies. To select a small sample of isolated massive 
galaxies in the NS A catalog with AA > G12 


created a modified Tinker et al. (2011) group catalog 
where the central galaxies m the group catalog were de¬ 
fined to be located at least three times their halo radius 
from another central galaxy. We define our high mass 
isolated criteria by successfully reco vering these isolated 
massive galaxies from the modihed Tinker et al. (2011| 
group catalog. 

Eor high mass galaxies with M* > 10^-^Mq, we sim¬ 
ilarly measure the 2P projected separation to the near¬ 
est more massive galaxy by at least 0.5 dex in stellar 
mass within a 2P projected distance of 7 Mpc and within 
1000 km s“^ in redshift. We call this distance dhost.o. 5 - 
We place the same cut on dhost,o .5 of 1.5 Mpc as the low 
mass sample. 

The cut on dhost,o .5 is not applicable to the most mas¬ 
sive galaxies in our sample because too few or no mas¬ 
sive galaxies exist above 0.5 dex in mass. Therefore, we 
use the SPSS and 2MASS catalogs to measure the fifth 
nearest neighbor surface density (E„ = n/irr'^), where n 
is smaller if less than five neighbors are found within 7 

















































































3 


Mpc. To be isolated, we require S„ < 1 Mpc~^. This 
cut on is motivated by the modified Tinker et al.| 
(20111 group catalog, since we maximize the recovery of 
this catalog by cutting on < 1 Mpc“^. Decreasing 
the cutoff in E„ does not appreciably effect our results. 
Removing this cut entirely will affect our results by in¬ 
cluding one gas-depleted galaxy at h igh masses relative 
to the rest of the sample (see § 3.21. Therefore to be 
considered isolated, we require a massive galaxy to have 
both dhost. 0.5 >1.5 Mpc and E„ < 1 Mpc“^. 

Our isolated low mass galaxies may have a nearby 
neighbor with simil ar stellar mass that ha s not been cap¬ 
tured by dhost (e.g. iStierwalt et al.||2015|). We calculate 
the 2D projected separation to the nearest, more mas¬ 
sive galaxy (dpair) for ah galaxies with M* < 10^° Mq. 
We use this environment variable to avoid radio beam 
confusion. We only require that dpair > 2rbeam where 
the beam radius (rbeam) is defined in Equation 10 be¬ 
low. Both halves of the overlapping pair are removed. 
We discard 53 galaxies, both isolated and non-isolated, 
with overlapping beams, all of which are from the AL¬ 
FALFA survey. We find no motivation to place further 
restrictions on dpair. In summary: 


(a) A low mass galaxy is isolated if M* < 10®'® Mq and 
dhost > 1-5 Mpc. 

(b) A high mass galaxy is isolated if M* > 10®'® Mq, 
dhost,o.5 > 1-5 Mpc and E„ < 1 Mpc“®. 


We refer to all galaxies that do not meet our isolation 
criteria as “non-isolated” for the duration of this paper. 

In terms of central and satellite galaxies: All isolated 
galaxies are central galaxies. All satellite galaxies are 
categorized as non-isolated. Due to the strictness of our 
isolation criteria, many non-isolated galaxies would be 
categorized as central galaxies. We define 36592 galaxies 
as isolated and 69782 galaxies as non-isolated in our cat¬ 
alog of 106374 total galaxies. We generally assume that 
our isolated sample is largely comprised of disk galaxies, 
even at the lowest masses. 


2.3. HI Observations 

We have obtained 317 HI observations with S/N 
greater than 2 in an observing program spanning 2005 
July to 2014 Sept, see G06 for an introduction to these 
data. After 2008, we focus our observations on iso¬ 
lated, low mass, highly inclined (an axis ratio less than 
0.65) galaxies from the NSA catalog described above with 
dhost > 1.5 Mpc, M* < 10®'®Mq. 

All single-dish, unresolved HI measurements were ob¬ 
tained using the 305 m Arecibo Telescope (AO) and the 
100 m Green Bank Telescope (GBT). For observations at 
AO in 2005 (hereafter A02005), we observed with the L- 
Band Wide receiver using 1024 resolution channels with 
a bandwidth of 12.5 MHz and a velocity resolution of 
2.6 km s“^ . For observations at AO in Spring 2013, 
Spring 2014 and Fall 2014 (hereafter AO2013, AO2014a 
and AO2014b), we observed using the L-Band Wide re¬ 
ceiver using 4096 resolution channels with a bandwidth 
of 12.5 MHz and a velocity resolution of 0.65 km s“^ . 
For observations at GBT in 2005, 2006 and 2008 (we 
group all observations together as GBT2008), we ob¬ 
served using the L-band receiver using 8192 channels 


with a bandwidth of 12.5 MHz and a velocity resolu¬ 
tion of 0.32 km s“^ . In all cases, position switched 
ON/OFF observing and ON/OFF noise-diode calibra¬ 
tion was used. 

Our target galaxies have an average Petrosian 90% 
light radius (rgo) of ~ 20", significantly smaller than 
the GBT and AO beam sizes (9' and 3.5', respectively), 
therefore flux attenuation and pointing offsets should not 
significantly affect our observations. We remove galaxies 
from our sample in the rare instances where the beam 
size is less than rgg. We detail each observing run in 
Tableincluding average integration times, the number 
of reduced spectra, the number of high S/N spectra, the 
number of non-detections and the number of unstable 
spectral baselines. For 371 successfully reduced spec¬ 
tra, we make 317 detections of galactic HI emission. We 
do not detect HI emission in 54 cases, most of which 
are non-isolated galaxies fr om G BT2008. We discuss our 
non-detections further in § |2.7[ 

2.4. HI Data Reduction and Measurements 


While there exist many methods for determining the 
HI properties of 21-cm emission lines, we implement a 
simple and consistent method optimized to minimize user 
intervention and maximize reproducibility of our mea¬ 
surements. We measure the shape and total emission of 
each single-dish HI observation without assuming a func¬ 
tional form for the 21-cm line. We also automate our 
entire procedure following a preliminary quality-control 
check. Automating our measurement procedure enables 
us to use a Monte Garlo method for e stimating random 
errors in our HI parameters (see § 2.51. 

We first use our observing notes to identify 23 spectra 
with bad spectral baselines which can significantly affect 
our HI measurements. If the signal cannot be differen¬ 
tiated from a noisy baseline then we discard the obser¬ 
vation and add the galaxy to a future observing queue. 
Next, we perform a linear baseline fit to each polarization 
mode outside of the region of each galaxy’s HI emission 
at the SDSS systemic velocity. Since we generally target 
isolated galaxies, we do not find overlapping signals in 
our spectra. We subtract the baseline fits from each po¬ 
larization mode and average the two polarization modes. 

To achieve higher S/N and ensure reliable detections 
with our measurement algorithm, we determine an opti¬ 
mal smoothing length for each observing resolution that 
produces the most accurate HI parameters in our reduc¬ 
tion pipeline. We test for optimal smoothing lengths by 
degrading synthetic Gaussian signals over a range of real¬ 
istic widths, integrated fluxes and signal-to-noise levels. 
We use Hanning windows to smooth spectra to a resolu¬ 
tion of 5 km s“^ for all of our observations. 

For each spectrum, we identify the systemic velocity 
channel using the optical SDSS redshift velocity. If a 
spectrum does not have a signal-to-noise greater than 
2 within 100 km s“^ of the systemic velocity, we flag 
the observation as a non-detection. We discuss non- 
dete ctio ns as well as the sensitivity of our observations 
in§[^ 

Next we identify the peak flux, /peaki within a 
500 km s“^ spectral window around the systemic ve¬ 
locity. We expand outward from the peak flux velocity 
channel and set the boundaries of the emission to the 
first locations where the spectrum drops below zero flux. 
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We locate these boundaries using a copy of the spectrum 
that has been smoothed by twice the Hanning window 
smoothing length. This extra smoothing eliminates any 
spurious noise features within the HI emission that may 
drop below zero and cut the emission profile artificially. 

We measure the velocity width at 0.5/peak, bbsoi by 
measuring the distance in velocity space between the first 
and last spectral channels with flux above 0.5/peak- To 
ensure we have measured the most accurate value of W 50 , 
we interpolate between the neighboring spectral channels 
on either side of the channel closest to 0.5/peak- 

We measure W 50 in order to directly compare with 
published ALFALFA results. However, since we have 
no reason to assume that W 50 is the best probe of the 
maximum circular velocity of all disk galaxies, we also 
measure W 20 , the velocity width at 0.2/peak- While W 20 
is more sensitive to noi se in the wings of the HI profile 
(Koribalski et al. 2004), our high S/N isolated galaxies 
may otter us a deeper dynamical probe into the HI disks 
than W 50 and VF 20 - We therefore focus on W 20 in our 
analysis below. We note that H 20 is typically 25 km s~^ 
larger than Hro h i our isolated sample, similar to |Korib-| 


Equation 2 of Haynes et al. (2011), 


alski et al. (|2004'1). 

We compute the integrated flux density, 5 ' 2 i, by inte¬ 
grating the flux between the channels of the emission line 
profile boundaries. In order to remove excess flux due to 
noise in the wings of the profiles, we fit a second order 
polynomial to the minima on either side of the emission 
profile between 0.5/peak an d where the flux goes to zero, 
similar to the algorithm of Springob et al. (2005). Using 
these fits, we remove spurious signal trorn the wings of 
the profile that tends to positively bias our flux measure¬ 
ments. We have tested this method by degrading high 
S/N spectra and then measuring the flux and HI pro¬ 
file widths as the S/N decreases. We found this method 
effective down to S/N of 2. 

We calculate the root mean square (rms) of the noise 
of each spectrum, cTrms, by identifying the largest signal- 
and RFI-free spectral chunks on either side of the 
HI emission profile with minimum window widths of 
200 km s“^ . Two windows are selected on either side 
of the emission line, beginning zLW^o /2 away from the 
profile edges and extending out to the maximum extent 
of the spectral range. Each window range is then opti¬ 
mized by shrinking the window until the mean flux is less 
than 0.5 mJy or the minimum window width is achieved. 
The restriction on the mean flux ensures that we select 
a chunk of the spectrum where the baseline is stable and 
no significant emission is contaminating our noise mea¬ 
surement. We track the mean flux in the signal-free and 
REI-free portion of the spectrum to monitor baseline sta¬ 
bility. We eliminate REI by masking any single channels 
with 6 -sigma emission or greater. We calculate the rms 
as the root mean square of the two signal-free, REI-free 
and baseline stable spectral chunks. 

If multiple observations exist for a galaxy at the same 
telescope and observing run, each rms weighted spectrum 
is co-added together. Eor multiple observations of the 
same galaxy over different observing runs, we compute 
the error weighted mean of each HI parameter measured 
with our reduction method. We tune and test our HI 
measurement method with the ALEALEA40 data in S 


S/N = 


W; 


50 


1/2 

U^smo 


( 1 ) 


with each variable as defined above and where 
'Wsmo = kF 5 o/ 10 [km/s]. Hanning smoothing redis¬ 

tributes spectral noise and therefore increases the signal- 
to-noise ratio of our spectra. This smoothing creates 
dependancies in each channel on neighbo ring channels. 
This d ependance is accounted for in Wsmo- [Haynes et al.| 
( 2011 ) calculates Wsmo as half of the number channels 
across the smoothed 50% HI line width. So iCgmo rep¬ 
resents the number of independent channels in each HI 
line. 

2.5. HI Parameter Error Estimation 

To estimate random errors in our HI parameter mea¬ 
surements, we apply a Monte Carlo (MC) technique. Eor 
each MC iteration, we add a random noise realization to 
the fiducial observed signal using our estimate of the rms 
noise in each spectrum and then remeasure the HI pa¬ 
rameters over 1000 trials. The errors of each HI param¬ 
eter estimate are the rms of the difference between each 
MC realization and the fiducial value of the parameter 
estimate. 

We note that our method of error estimation does not 
consider baseline stability and its influence on the inte¬ 
grated flux density, which may be a source of significant 
uncertainty. It is for this reason that we have visually 
inspected and removed any spectra with obvious base¬ 
line instability. We have also made an effort to obtain 
multiple observations of the same galaxies to inspect any 
add itional systematic error (see the last paragraph of § 


2.6 We calculate a S/N measurement for all of our ob¬ 


servations following the ALFALEA survey, adopted from 


2.6. Homogeneous HI Measurements of ALFALEA 

In addition to our observations of isolated low mass 
galaxies, we have re-measured all Arecibo Legacy East 
ALFA spectra overlappi ng with our NSA ca talog (from 
the 40% data release of Haynes et al. (2011)). The AL¬ 
FALFA survey is a blind, single-dish, hux-limited, HI sur¬ 
vey. The ALFALFA40 data release covers ~2800 deg^ of 
the sky and provides a significant amount of overlap with 
the SDSS footprint. We match ALFALFA HI sources 
where the HI center is within 25" of the optical center 
of our NSA catalog, given the pointing accuracy of the 
survey ( Haynes et al.jpMT ). We use the ALFALFA data 
to compliment our isolated low mass sample with both 
non-isolated galaxies and massive galaxies. 

Re-measurements are preformed primarily to maintain 
a homogeneously measured data set and to apply our 
MC error calculation. We do not visually inspect or 
baseline fit the ALFALFA40 data, instead we begin with 
the baseline-corrected spectra from ALFALFA and cut 
on their quality code of 1 to ensure reliable detections. 
We run these ALFALFA40 data throu gh o ur HI param¬ 
eter measurement steps described in § 2.4 convolving a 
3 channel Hanning smoothing window with the down¬ 
loaded ALFALFA40 spectra. 

Figure shows the difference between our re¬ 
measurements and the original ALFALFA40 measure¬ 
ments in units of the error of each HI parameter between 
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Fig. 1. — Two dimensional distribution of the difference between 
our re-measurement of the ALFALFA spectra and the original AL¬ 
FALFA measurement in units of a, where cr is the error for each 
galaxy from the two data sets added in quadrature. The difference 
in 50% velocity width Wso and HI flux density S 21 are plotted 
on the X- and y-axis respectively. The number of galaxies in each 
bin are shown in the greyscale bar on the right. The over plotted 
box represents our selection boundaries (di2cr) on the data. Any 
measurement outside of the box is discarded as a failed measure¬ 
ment. We discuss the asymmetry in the distribution at negative 
differences in W 50 and positive differences in 521 which is caused 
by a difference in measurement methods in § |2.6| 


the two data sets added in quadrature. In general, both 
our remeasurements of W 50 and 521 agree extremely well 
with the ALFALFA40 measurements. As expected, we 
find that just under 5% of galaxies lie beyond the ±2cr 
level (4098 HI observations in Figure!^ (169 outliers lie 
beyond the 2.0-cr boundaries in either IA50 or 52i). We 
remove galaxies beyond the 2 (t window from our sam¬ 
ple because their emission line boundaries are ambiguous 
and we require minimal scatter in our analysis below. We 
also discard any galaxy where the error on either IF50 or 
Mhi is two times the measured value or more. We re¬ 
cover 96% of the complete ALFALFA40 dataset with our 
remeasurements. 

We note an asymmetry in the 2D distribution of differ¬ 
ences in W 50 and 52i in Figurel^where we underestimate 
W50 and overestimate 52i relative to ALFALFA. The 
cause of this is due to a difference in measurement meth¬ 
ods. Since we measure IF50 at 50% the peak flux with¬ 
out accounting for noise or the shape of the HI profile, 
our value of W 50 will always be larger than ALFALFA 
since they account for the shape of the HI profile when 
calculating such line-widths. This will tend to decrease 
the value of the peak flux and therefore ALFALFA will 
always measure an equal or larger value of W^o- The 
spectral boundaries for the A LFALFA measuremen t are 
manually marked by eye (see Haynes et al. (2011), § 3 
item 6). It is possible that this tends to remove some 
of the actual signal that overlaps with the noise features 
in the wings of the HI profile while our method of auto¬ 
matically marking spectral boundaries tends to include 
some fraction of the noise in the wings of our profiles. 
Since less than 1% of our measurements are beyond ±lcr 
in this region, and 3% are beyond ±0.5(T, we proceed 
with our measurements. We also note that our analysis 
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Fig. 2.— A sample pair of overlapping HI spectra between our 
sample and ALFALFA40. The left spectra are from our observa¬ 
tions. The right spectra are from the ALFALFA40 data set. In all 
panels, the smoothed spectrum is plotted as a grey dotted line, the 
emission region is identified in solid black, the 50% and 20% flux 
boundaries are identified with red dash-dot and blue dash-dot-dot 
vertical lines respectively. The central velocity is marked with a 
vertical dotted black line. The top spectra are representative of a 
straight-forward measurement, while the bottom spectra are more 
ambiguous due to the difference in spectral resolution. The more 
ambiguous cases will result in larger errors in W 50 and 521 due to 
our method of estimating errors. 


proceeds with IF20, which is less sensitive to noise in the 
peak flux and a more effective measurement of line-width 
in the lowest mass galaxies. We will discuss this further 
in our future paper on the baryonic Tully-Fisher (BTF) 
relation. 

Our HI observations overlap with 15 of the AL- 
FALFA40 sample. We compare our HI parameter esti¬ 
mates to the original ALFALFA40 catalog for each over¬ 
lapping galaxy. Most of our HI parameter estimates 
agree when compared to the original values from the AL- 
FALFA40 data release. However, several measurements 
of S 21 disagree between our GBT2008 observations and 
the original ALFALFA40 catalog. We assume that both 
our measurements and the ALFALFA40 measurements 
are correct and we attribute any differences to system¬ 
atic errors that we have underestimated with our error 
estimation method (most likely due to baseline instability 
in the GBT2008 data). In order to bring our GBT2008 
measurements into agreement with the ALFALFA cat¬ 
alog, we add 0.2 Jy km/s in quadrature to all of the 
GBT2008 S'21 uncertainties. We present two pairs of 
overlapping galaxy spectra in Figure for reference. In 
this figure, the top pair is a straightforward comparison 
while the bottom pair is more ambiguous. 

2.7. Non-Detections 

Our observing strategy using the AO is to integrate 
until we measure a significant HI signal or until we reach 
a 25 minute on-source integration time when possible. 
Most of our observations result in detection, but in some 
cases we do not detect significant HI emission. We de¬ 
fine a non-detection as an observation with S/N < 2, 
the threshold where our measurement method begins to 
fail. Out of our 371 reduced spectra, 54 have a S/N less 
than 2. Nearly all of these low S/N observations are from 
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GBT2008 since these observations were made before we 
implemented a fixed observing time criteria and focused 
our observations on isolated galaxies. Given the smaller 
dish size, the average GBT exposure leads to smaller S/N 
than an equivalent AO integration. Out of 144 isolated 
galaxies we observed using our post-2008 observing strat¬ 
egy, 1 are non-detections. 

For all HI non-detections, we determine upper limits 
on 5'2 i. First, we make preliminary estimates of atomic 
gas mas ses using the stellar-to-atomic gas mass relation 
in § 3.2 Next, we estimate the inclination corrected 50% 
velocity widths using the relation between baryonic mass 
and 50% inclination-corrected veloc ity widths, similar to 
the relation we describe in § 3.3.1| We inclination de- 


correct this 50% velocity width to obtain an estimate 
of Wso- We initialize the value of 521 assuming a gas 
fraction of /gas = 0.98. 

Next, we create a synthetic Gaussian 21-cm emission 
line with a FWHM equal to W 50 and a total integrated 
flux of 521. We inject this initial Gaussian signal into 
100 noise realizations using the measured Urms of our 
observation. After each noise realization, we reduce 52i 
by 10% and repeat a new set of noise realizations while 
we keep the FWHM of the Gaussian fixed. We repeat 
this procedure until the S/N drops below 2. The final 
injected values of 52i are reported as the upper limits 
of each non-detection. These upper limit estimates are 
plotted as light-colored, downward pointing arrows in all 
plots. Upper limits are excluded from our velocity-mass- 
size scaling relations and baryon fraction analysis since 
their velocity widths are contrived. 

ALFALFA40 code 1 sources are all above the AL¬ 
FALFA survey’s S/N threshold of 6.5, therefore upper 
limit calculations are unnecessary. Gode 2 sources from 
ALFALFA40 have lower S/N and since they do not sig¬ 
nificantly impact our results below, we do not consider 
these data in our analysis. 

2.8. Calculated and Inferred Quantities 

We next describe the quantities derived directly from 
our HI line emission parameters and the optical param¬ 
eters from the NSA catalog. In Table [l] we present a 
small sample of isolated low mass galaxies from both 
our 255 original observations and the re-reduction of 
4032 ALFALFA40 observations. Our complete data set 
is available upon request and will be publicly released 
once our HI observation program has concluded. 

Distances {D) reported here are luminosity distances 
calculated from pe culiar motion corrected redshifts 
(Willick et al. 1997). For our calculations and analysis, 
we assume that all of our galaxies are oblate, axisym- 
metric spheroids with symmetric mass distributions. We 
infer the observed inclination, sin i, from the axis ratio 


as, 


sm i = 


'1 - [hjaY 


( 2 ) 


where is the intrinsic axis ratio, b/a is the observed 
axis ratio and i is the inclination angle. Note that we 
limit the maximum value of the axis ratio b/a to 0.9995 
to avoid dividing by zero in Equation!^ 

The value of is traditionally fixed at about 0.2 at 


all stellar masses ( e.g., G06, Begum et al. 2008). See 
Yuan & Zhu| (2004) for a brief discussion on the history 
of go as well as the distribution of go as a function of 
morphology. The effect of go on inclination corrections 
can be significant enough to change the slope and pos¬ 
sibly affect the scatter of the BTF relation, (see Singhal 
2008| Figure 1.8). We plan to discuss the effect of go 
further in a future paper, where we will analyze circular 
velocities using the traditional value of go = 0.2 as well 
as with a prescription for a changing go as a function of 
stellar mass. In this work, we set go = 0.2 throughout 
our analysis below. 

We remove redshift broadening and we inclination cor¬ 
rect the HI line-widths to estimate a measure of maxi¬ 
mum rotational velocity of the stars and HI gas disk in 
each galaxy as. 




Wx 


2sin i(l -I- z) ’ 


( 3 ) 


with X = 50 or X = 20 for either the 50% or 20% 21-cm 
line-width, respectively. We perform these inclination 
corrections to Wx but we do not correct for non-circular 
motions (i.e. turbulence) which is typically assumed to 
be small (less than 20 km s“^ ). Note that in previ¬ 
ous work, G06 corrected for non-circular motions. We 
note that the 20% velocity width has been shown by 
some authors to provide a better estimate of t he true 
maximum circular velocity of the galactic disk (Verhei- 
jen 1997 Avila-Reese et al. 2008). We plan to discuss 
the effect of non-circufar motions on HI velocity widths, 
as well as the accuracy of HI line-widths as a probe of 
maximum circular velocity in a future paper. For now, 
we assume Vw 20 ,i is an accurate representation of maxi¬ 
mum circular velocity. 

We find that the inclination measurement for a small 
fraction of galaxies in our catalog fail due to the pres¬ 
ence of central galactic bars. Inclination measurements 
fail when a face-on galaxy resembles the light distribution 
of a more edge-on galaxy due to the thin, centrally con¬ 
centrated distribution of light of a galactic bar. We re¬ 
move any barred gal axies that have bee n identified by the 
Galaxy Zoo project (Hoyle et al. ]MT| ) from our analysis 
of circular velocities. We preserve these barred galaxies 
unless specified in our analysis below. 

We calculate the HI mass from the 521 flux us ing the 
standard formula (e.g. Haynes fc Giovanelli||1984) as, 


M„, = 2.356 X 10= (X) (4) 

which assumes the HI gas disk is optically thin. We do 
not correct for HI self-absorption so our 21-cm flux mea¬ 
surements are lower limits due to the self-absorption of 
HI flux in a galaxy (Haynes & Giovanelli 1984 Haynes 

eFalllMI] ). 


We present HI masses as a function of redshift in Fig¬ 
ure This figure shows the ALFALFA40 dataset as grey 
points with our data and data from GOG plotted in pur¬ 
ple. Our observations probe to lower HI masses at all 
redshifts than ALFALFA40. 

We calculate atomic gas masses as Mgas = 1.4Mhi 
where the factor of 1.4 corrects for helium assuming so¬ 
lar abundances. We calculate “cold” baryonic masses 
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Fig. 3. — HI masses versus peculiar motion corrected red- 
shifts/distances of NASA Sloan Atlas galaxies that have been ob¬ 
served in 21-cm emission. Grey data points are our ALFALFA 
re-measurements. Purple squares are our own new and previ¬ 
ously published (G06) observations with downward pointing, light- 
colored arrows representing upper limits on HI masses for non¬ 
detections. Our HI observations complement the ALFALFA sur¬ 
vey at low masses and also overlap the ALFALFA sample in several 
cases as a consistency check between datasets. 

as Af^baryon — ^ga.s 

fraction defined as. 


M* with the atomic gas mass 


/gas — 




gas 


Ml 


baryon 


(5) 


Note that our final cold baryon masses do not include 
dust, molecular gas or ionized gas, the combination of 
which may cont ribute significantly to the baryon content 
of our galaxies (IMc Gaugh fc de Blok[|1997| ICortese et al. 
2012 Gnedin|20l2p. f urther quantities, mcludmg stellar 
mass surface density and baryon fraction are described 
in the appropriate sections below. We discu ss th e impact 
of molecular gas on our results in section § 4.1 


2.9. HI Detection Limits, Isolated vs Non-Isolated 
Sample Sizes and Surface Brightness Limits 

Here we discuss the relative sample sizes of our isolated 
and non-isolated galaxies as well as the surface brightness 
limitations of the NSA catalog. We also discuss how a 
flux-limited HI survey affects our results and subsequent 
conclusions. 

2.9.1. HI Detection Limits 

The minimum HI mass detected by ALFALFA that 
overlaps with our catalog is 10 ^'^Mq while the minimum 
HI mass detected by our observations is lO^’^M©. Since 
the ALFALFA survey is flux-limited, it has not detected 
low mass, Hl-depleted galaxies such as satellites in dense 
groups and clusters (see § 6 of Haynes et al. (2011) for a 
thorough discussion of the ALlALlAdU completeness). 
For the duration of this paper, we define detection limits 
as the minimum HI mass of either our data set or the 
ALFALFA40 catalog that overlaps with the NSA catalog. 

When the ALFALFA survey is combined with our ob¬ 
servations, the HI sample we study here is largely incom- 
plet e and biased towards blue, Hl-rich, lat e-type galax¬ 
ies (Haynes et al.|2011 |Huang et al.||2012a).The median 
q — r color tor the AlIaLIA survey and our HI data 


is 0.5 and 0.3, respectively. Sat ellite galaxies in the lo¬ 
cal group dSpekkens et al.||2014|) and in the Virgo clus¬ 
ter ( Cortese et ai.||2011 ) are observed to be severely HI- 
depleted and would not be detected by either our obser¬ 
vations or by ALFALFA. Therefore, even the relatively 
Hl-depleted galaxies in our sample are still quite Hl-rich 
compared to the average red sequence galaxy. 

Due to our HI observation strategy, we have detected 
nearly every isolated galaxies identified optically. For the 
few isolated galaxies that we have not made a detection, 
we calculate upper limits at relatively large atomic gas 
fractions. This is one of the primary results of our paper: 
all of the isolated galaxies that we have observed have ei¬ 
ther been measured at large atomic gas fractions or large 
atomic gas fract ions cannot be ruled out via their upper 
limits (see §3.1). As for our comparison to non-isolated 
galaxies, our sample only glimpses into the Hl-poor pop¬ 
ulation that we expect to find in satellite galaxies. In our 
analysis below, we discuss these points further in context 
with our results. 

2.9.2. Isolated versus non-isolated galaxy sample sizes 

In order to disentangle internal galaxy evolutionary 
processes from the effects of environment, we select only 
the most isola ted g alaxies in the NSA catalog given our 
definition in § 2.2 Despite our overall HI sample being 


largely incomplete, we have made every effort to observe 
similar numbers of isolated and non-isolated galaxies at 
low stellar masses. Our analysis therefore assumes that 
the outliers we measure are significant and not simply 
an effect of the relative sample sizes between isolated 
and non-isolated galaxies. 

For galaxies with stellar mass below we 

have 429 non-isolated galaxies and 368 isolated galax¬ 
ies with HI observations. We achieve an average frac¬ 
tion of isolated galaxies with HI observations of 0.46 for 
M* < while the average fraction of isolated 

galaxies in the entire NSA catalog is 0.34. We assume 
that this relative sample size should allow us to rigor¬ 
ously compare the relations, distributions and outliers as 
a function of environment. 

2.9.3. Surface brightness limitations 

Even though the NSA catalog has been optimized for 
nearby galaxies in the SDSS, the catalog still suffers from 
low surface brightness completeness limitations. At a 
half-light surface brightness /i5o,r ~ 23.5 mag arcsec”^, 
the SDSS completeness drops below 50%. At ^50,r ~ 
24.0 mag arcsec~^ the SPSS com pleteness drops below 
10% dBlanton et al.||2005b| [20081 G12). 

This surface brightness incompleteness impacts our 
study if either of two extreme situations exist: a.) we 
preferentially miss Hl-poor galaxies or b.) we preferen¬ 
tially miss gas-rich galaxies due to incompleteness. 

If we preferentially miss Hl-poor galaxies, then a Hl- 
poor population of isolated galaxies might exist that has 
gone undetected. G12 compared the surface brightness 
distributions of low mass non-isolated galaxies to isolated 
galaxies via the two-sided K-S test and found that the 
two samples are likely drawn from the same distribution. 
Therefore, we have little evidence that suggests a larger 
number of Hl-poor isolated galaxies at the low surface 
brightness end of the distribution are missing from our 





































sample compared to the number of Hl-poor non-isolated 
galaxies. 


The recent discovery of ultra diffuse galaxies (UDGs) 

in the Goma cluster by 1 

yan Dokkum et al. 

2015a|b|) 

using the Dragonhy Telephoto array 1 

Abraham & van 


icant population of galaxies with extremely low central 
surface brightness exists in galaxy clus ters and the local 
field (24 < < 26 mag arcsec”^) (Dalcanton et al. 


19971. This may indicate that more low mass galaxies 
(with —13 < My < —17) are missing from our SDSS 
sample than previously thought. It is possible that a sim¬ 
ilar population of UDGs exist in “the held”. However, we 
have no evidence at this time that we are preferentially 
missing a significant number of isolated Hl-poor UDGs. 

We have no evidence that a signihcant number of HI- 
rich LSB galaxies are missing from our sample due to the 
following logic: We hnd that galaxies with surface bright¬ 
nesses below fir ,50 < 2 3 mag arcsec~^ have a m edian /gas 
of 0.7 ± 0.2 (also see Schombert et al. (2001)). There¬ 
fore, at all stellar masses where fir ,50 < 23 mag arcsec“^, 
galaxies should have atomic gas masses equal to or signif¬ 
icantly greater than their stellar masses. These sources 
would be detected by ALFALFA but would have no op¬ 
tical counterparts in SDSS. The ALFALFA40 data re¬ 
lease has found just ~ 50 extragalactic HI sources (out 
of 15,855) without optic al counterparts in the SDSS DR7 
and a ssociated catalogs (Haynes et al. 2011 Cannon et al. 


2015). This suggests that while such a population exists, 
these “dark galaxies” are rare. 

3. RESULTS 

We begin our analysis by inspecting the atomic gas 
fraction of galaxies as a function stellar mass and stellar 
mass surface density. We then compare the relation and 
scatter of stellar mass to atomic gas mass in isolated envi¬ 
ronments compared to non-isolated galaxies. We analyze 
the velocity-mass-size galaxy relations and the impact of 
environment on the slope, scatter and correlation of these 
scaling relations. Finally, we investigate baryon fractions 
in these systems. 

In Figures 1^ through we represent isolated galaxies 
with red squares and non-isolated galaxies with black 
dots. Whenever we exclude galaxies from our analy¬ 
sis, we preserve these data as grey points. We plot up¬ 
per limits as light-colored, downward pointing arrows, 
color coded with environment. We also plot a hollow 
data point with a representative set of error bars wher¬ 
ever error bars are not explicitly presented. In Figures 

and [7 we plot non-isolated galaxies we determine to 
oe gas^epleted as open blue dots. We hav e applied the 
Bayesian linear regression fitting method of Kelly (20071 
for all linear fits below. This fitting method is particu- 
larly useful for linear models that also measure the in¬ 
trinsic scatter of galaxy scaling relations. 

3.1. Atomic Gas Fraction Versus Stellar Mass 

The fraction of baryons in atomic gas, /gas, is a st rong 
indicator of the evolution ary state of a galaxy (e.g., Mc- 


Gaugh & de Blok 1997). Spheroidal, low mass galax¬ 
ies m dense environments can be almost completely de¬ 
pleted of HI gas, s uch as the satellite ga laxies of the 
Local Group (e.g., Spekkens et al. 2014). Low mass 
galaxies in moderately isolated environments have been 


shown t o form stars slowly and in sparsely distributed 
regions (van Zee|2000 2001) with star formation bec om- 
ing more stochastic at low inasses ( Kauffmann|2014 ). In 
fact, the timescale for gas consumption due to star for¬ 
mation alo ne in gas-rich , irregular, low mass galaxi es is 
^ 20 Gyr (|van Zee||2001| |Roychowdhury et al.||2014D . 

Given the timescale tor gas depletion due to star for¬ 
mation in low mass galaxies and the fact that isolated 
low mass galaxies are all forming stars (G12), isolated 
galaxies should contain a significant amount of baryons 
in atomic gas. We might also expect gas-depleted low 
mass galaxies to be located in dense regions and to have 
undergone environmental processes, as opposed to have 
experienced purely internal processes that might deplete 
atomic gas. Therefore, examining the atomic gas frac¬ 
tions of galaxies in isolation compared to denser envi¬ 
ronments can provide us with insight into how feedback 
processes and environmental effects regulate the atomic 
gas supply of galaxies. 

We present the relation between stellar mass and 
atomic gas fraction in the left panel of Figure]^ and the 
relation between stellar mass surface density and atomic 
gas fraction in the right panel of Figure]^ 

In Figure0(Left), at stellar masses above ~ IO^'^^Mq, 
galaxies follow a scattered relation with an upper ridge- 
line on /gas at all stellar masses. At stellar masses below 
109.25^0, atomic gas fraction appears to saturate 
at M* ~ 10®’^® with /gas ~ 95%. This implies that 
high mass galaxies are either more efficient at converting 
their gas into stars or are more effective at preventing 
gas from cooling than galaxies at low masses. The scale 
dependency on /gas at low masses is not introduced sim¬ 
ply because of our dynamic range in this figure. This 
fact is made more apparent by Figure ^ where a break 
in the relation between stellar mass and atomic gas mass 
is clearly observed. 

While the median atomic gas fractions for isolated and 
non-isolated low mass galaxies are similar (0.81 ±0.13 for 
isolated galaxies and 0.82 ± 0.16 for non-isolated galax¬ 
ies) isolated galaxies present slightly less scatter in /gas 
at fixed stellar mass. Indeed, non-isolated galaxies with 
M* < 10 ®'^®Mq here and satellite g alaxies from the 
litera ture (e.g.. Galactic dSphs from Spekkens et al. 
2014), can inhabit nearly the full range of /gas (from 


U.OU to 0.95) while we observe all isolated galaxies with 
M* < 10 ®-^®Mq to have /gas > 0.3. 

We test that the lack of isolated galaxies with /gas < 
0.3 is not simply due to a sampling effect. Over 10,000 
trials, we randomly select the number of isolated low 
mass galaxies (A^isoiated = 368) with replacement from all 
797 galaxies with M* < 10® On average, we draw 

5 galaxies with /gas < 0.3 and we always draw at least 
1 galaxy with /gas <0.3. The complete lack of observed 
isolated galaxies with /gas < 0.3 is therefore compelling 
evidence for the effects of environmental processes on 
atomic gas regulation. 

We observe atomic gas fractions in our low mass iso¬ 
lated sample that vary from 0.30 to 0.99. Therefore iso¬ 
lated galaxies retain a significant fraction of baryons in 
atomic gas at low masses. The fact that we do not de¬ 
tect isolated low mass galaxies with /gas below 0.3 is a 
strong constraint on internal feedback processes, since it 
appears that feedback processes alone in low mass galax- 
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Stellar Mass vs Gas Fraction Stellar Surface Density vs Gas Fraction 
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Fi g. 4 . — (Left) Atomic gas fraction versus Stellar mass. Black dots are galaxies that hav e no t passed the isolation criteria described 
in § |2.2| Red squares have passed the isolation criteria for each mass regime described in § |2.2| We plot the median error bar for the 
entire sample as the open point, which is representative of the rest of the sample. This relation saturates at about 10^ Mq in stellar 
mass. The dashed grey line represents the detection limit (minimum HI mass detected) of the overlap of our sample with the ALFALFA 
survey (10 ^'^Mq). The dashed purple line represents the detection limit (minimum HI mass detected) of our HI observations (IO^'^Mq). 
Non-isolated galaxies are scattered to lower /gas compared to isolated galaxies at nearly all stellar masses, indicating that environment 
affects the atomic gas fraction of non-isolated galaxies. Low mass isolated galaxies are not found below an /gas of 0.3. (Right) Atomic gas 
fraction versus stellar mass surface density. (Right) We divide stellar mass by the square of the effective radius to show that stellar mass 
surface density is a better predictor of atomic gas fraction than stellar mass alone. 
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Fig. 5. — (Left) Atomic gas mass plotted as a function of stellar mass. Black dots are non-isolated galaxies, red squares are isolated 
galaxies. We plot the median error bar for the entire sample as an open point, which is representative of the rest of the sample. The fits 
here are listed below and in Table and are plotted as the solid black line. We plot a dashed blue line representing the —1 dex threshold 
where we define galaxies to be gas-^pleted when located below this threshold. We identify non-isolated gas-depleted galaxies as open blue 
dots. (Right) Residuals to our fit as a function of stellar mass for all non-isolated galaxies (Top, Right) and all isolated galaxies (Bottom, 
Right). The break in the power-law is shown as a vertical, dashed grey line. The rms scatter of the residuals is listed at the bottom left 
corner of each panel. We plot a dashed blue line representing the —1 dex threshold where we define galaxies to be gas-depleted when located 
below this threshold. Note that no isolated galaxies are found below a residual of —1 dex. We plot non-isolated gas-depleted galaxies as 
open blue dots. The non-isolated galaxies are slightly more scattered but the primary relation between stellar mass and atomic gas mass 
is the same for both the non-isolated and isolated galaxies. 

ies does not remove all of a galaxy’s atomic gas. We note that the isolated galaxies with upper limit calcula- 
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tions cannot be conclusively ruled out as having lower 
atomic gas fractions than this limit. Feedback processes 
are most likely important in driving the large observed 
scatter in /gas at low masses. 

All but one non-isolated galaxy with M* < 10® Mq and 
/gas < 0.3 are observed by the blind ALFALFA survey. 
For clarity, we have over plotted a dashed grey line in Fig¬ 
ure 0 (Left) which represents the smallest HI mass that 
overlaps with the NS A catalog and therefore the smallest 
/gas that the ALFALFA survey could reach. Due to sig¬ 
nificantly longer exposure times, our observing program 
can detect galaxies at even lower atomic gas fractions 
than this grey line. We also over plot a dashed purple 
line in Figure]^ (Left) which represents the smallest /gas 
that our survey could have detected (the minimum HI 
mass detected by our observations is lO^^M©). We ob¬ 
serve one non-isolated galaxy and t wo is olated galaxies 
below the ALFALFA limit. Also see § 2.9 for a discussion 
of the sensitivity of our detections. 

Galaxies with large /gas are either less efficient at con¬ 
verting their gas into stars or these galaxies are accret¬ 
ing enough gas to replenish their star-making fuel sup¬ 
ply. As expected fro m galaxies that follo w the Kennicutt- 
Schmidt relation (see Zhang et al.||2009 section 3. 2) and 


from previous obse rvations (e.g. Gatinella et al. 2010 
Cortese et al.|20lT| ), we find that /gas anti-correlates with 
the surtace density distribution of stars. In the right 
panel of Figure we present the relation between stel¬ 
lar mass surface density and atomic gas fraction. We 
calculate the stellar mass surface density as. 


/i* 


d-Trr® ’ 


( 6 ) 


where reS is the effective radius which encloses 50% of the 
optical light. The effect of environment on /gas is essen¬ 
tially removed from this relation, indicating that galaxies 
may evolve along this relation. The distribution of stel¬ 
lar mass density is clearly a better predictor of /gas than 
stellar mass alone. The galaxies with the highest atomic 
gas fractions tend to have the lowest surface density of 
stars. 


3.2. Gas Mass Versus Stellar Mass 

We present the stellar-to-atomic gas mass relation in 
the left panel of Figure The slope of this relation 
steepens just below a stellar mass of 10 ®Mq. This break 
in the relation is the cause of the saturation of /gas at 
low mass in the left panel of Figure [5 

We fit a broken power-law to the hi^ and the low mass 
isolated galaxies only. The fit breaks at M* = IO^'^Mq, 
where we find the two linear fits best intersect and min¬ 
imize the scatter in the fits. This is a slightly lower stel¬ 
lar mass than our definition of a low mass galaxy. We 
present our fit s bo th in this panel and in Table As 
discussed in § 2.9 we find no evidence that this oreak 
in the power-law is due to a selection effect: We have no 
evidence that we have missed a significant galaxy pop¬ 
ulation at low stellar mass and high atomic gas masses. 
This break in the power-law has also b een observed by 
other studies (e.g., Huang et al. 2012a|. The break in 
this relation implies that star formation is either much 
less efficient in low mass galaxies or the HI gas in low 


mass galaxies stays hot enough to prevent star forma¬ 
tion but not hot enough to ionize the HI. 

We inspect the residuals to our fits in the right panel of 
Figure The residual panels demonstrate that scatter 
towards relatively small atomic gas masses is significant 
for non-isolated g alax ies but not so for isolated galaxies. 
As discussed in § |2.9[ the non-isolated panel is certainly 
incomplete at low atomic gas fractions due to ALFALFA 
survey limits, yet we still observe many galaxies with 
relatively large negative residuals. In fact, we only detect 
non-isolated galaxies at residuals less than —I dex. All 
isolated galaxies have either been either been detected in 
HI or cannot be ruled out at high HI masses via their 
upper limits. Therefore, we do not detect any isolated 
galaxies below a residual of —I dex. 

Note that —I dex is more than —3cr away from our pri¬ 
mary relation (see the rms values in each panel of Figure 

S (Right)) and we plot a dashed blue line in both panels 
this figure that represents residuals of —I dex. Be¬ 
low this residual threshold, we define gas-depleted galax¬ 
ies as open blue dots in this and our remaining figures 
for reference. We choose —I dex as our residual thresh¬ 
old because we do not observe any galaxies above I dex 
in either sample and we do not observe any isolated 
galaxies below —I dex. These panels reinforce the idea 
that environmental effects preferentially remove atomic 
gas from non-isolated galaxies. The non-isolated out¬ 
liers with residuals below —I dex provide evidence that 
galaxies are depleted of atomic gas due to environmental 
processes. 


3.3. Maximum Circular Velocity, Galaxy Mass and 
Galaxy Size Relations 

Distinct dark matter halos and their respective disk 
galaxy pr operties are related via galaxy scaling rela - 
tions (e.g.|Courteau et al.|2007t|Avila-Reese et al.|2008|). 
These relations are critical tor studying models of galaxy 
formation and evolution (e.g. van den Bosch||2000 Dut- 
ton et al.|2007 ). We measure the velocity-mass-size scal¬ 
ing relations or our isolated disk galaxies using HI veloc¬ 
ity widths, baryonic masses, and optical radii. Here we 
focus on galaxy scaling relations in isolation and compare 
our results to galaxies in denser environments. 

In each panel of Figure we present the velocity- 
mass-size relations as well as a fit to each relation using 
only unbarred, isolated galaxies with axis ratios less than 
0.65. For comparison, we also plot the fit for non-isolated 
galaxies only as a yellow dashed line in each figure. As 
mentioned in § |2.8[ we have removed known barred galax¬ 
ies from our analysis because their inclination measure¬ 
ments tend to fail. In Table we present the results 
of fitting a linear model of logy = a log a; -I- /3 ± e to 
our data, where e represents the intrinsic, random scat¬ 
ter in y. We also provide a measurement of the Pearson 
product-moment correlation coefficient (Pearson r) for 
each galaxy scaling relation. We provide model fits and 
Pearson r measurements for I.) all galaxies, 2.) isolated 
galaxies and 3.) non-isolated galaxies. 

We inspect the velocity-mass-size scaling relations of 
isolated disk galaxies. For maximum circular velocity 
measurements, we find that at Vw 50 ,i < 100 km s“^ 
the inclination corrected 50% velocity width can signifi¬ 
cantly deviate from the galaxy scaling relations of more 
massive galaxies, therefore we use the inclination cor- 
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slope=0.277+/-0.004 

const=-0.672+/-0.041 

scatter=0.075+/-0.002 
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Fig. 6.— Velocity-mass-size scaling relations for Vw20,i^ ^baryon reff. Unbarred, isolated galaxies with axis ratios less than 0.65 
are plotted as red squares. Non-isolated galaxies with axis ratios less than 0.65 and barred isolated galaxies are plotted as grey dots for 
reference. Gas-depleted, non-isolated galaxies with axis ratios less than 0.65 are plotted as open blue dots. The slope, constant, scatter and 
Pearson correlation rank (r) are listed in the top left corner of each panel. An open data point with representative error bars are plotted in 
the bottom right corner of each panel. The fit of each relation is plotted as a solid black line with ±1 (t scatter plotted as dotted black lines. 
We also plot the relations measured with the non-isolated sample only as dashed yellow lines for reference. The fits are similar for isolated, 
non-isolated and both samples combined (see Tablej^. The relation between Vw 2 Q,i ^-nd appears to be affected by environment. 

alone. Finally, we study the sizes of our sample using 
effective radii, VeS- Other choices of Petrosian isophotal 
radii are similarly scattered. 

3.3.1. Velocity versus Baryonic Mass 

In Figure (Top, Left) we show the relation between 
baryonic mass and maximum circular velocity, the BTF. 
We present two slightly different hts to this relation in 


rected 20% velocity width {Vw 20 ,i)- For mass measure¬ 
ments, Mbaryon has been shown to minimize the sc atter 


of the BTF for low mass, gas-dominated galaxies ( Mc- 
Gaugh et al. 1 120001 McGaugh] |2()05[ [Stark et al.| 20l)9[ 
Mcdaugh 2012 r ^baryon should therefore provide the 
best mass measurement for studying galaxy scaling re¬ 
lations, as opposed to optical luminosity or stellar mass 
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Table lousing mass and velocity as the independent vari¬ 
able in each case. This relation is the most linearly corre¬ 
lated with a Pearson r of 0.916 and only slightly decreases 
to 0.831 for the non-isolated sample. The classical Tully- 
Fisher relation is between optical luminosity and circular 
velocity while the BTF relates the total baryonic mass of 
the disk to the gravitational pot ential of the disk (Tully 
& Fisher 1977 Mo et al. 20101. When baryonic mass 
is used instead of luminous mass, galax ies at all masses 
follow a tight relation at all velocities (McGaugh et al. 
2000|. This indicates that despite being gas-dommated, 
low mass, irregular galaxies still follow the same fun¬ 
damental relation as more massive, stellar mass domi¬ 


nated galaxies (McGaugh et al. 2009). Indeed, we ob¬ 
serve a very small scatter ol 0.075 tor the isolated sam¬ 
ple and 0.080 for the non-isolated sample. We do not 
observe a significant slope change due to environment. 
Gas-depleted non-isolated galaxies are not preferentially 
scattered above or below this relation. We plan to ex¬ 
amine the BTF relation in more detail in a future paper. 


3.3.2. Radius versus Baryonic Mass 

In Figure (Bottom, Left) we show the relation be¬ 
tween radius and baryonic mass. The relation is linearly 
correlated with a Pearson r of 0.880 which decreases to 
0.764 for non-isolated galaxies. This relation is an esti¬ 
mate of the overall density and thus the c oncentration 
of a galaxy. Firmani & Avila-Reese (2009) have shown 
that disk galaxies move along this relation as they evolve, 
indicating that we may not observe environmental ef¬ 
fects on the baryonic mass to radius relation if radius 
and baryonic mass are both modified by environmental 
processes. Indeed, while non-isolated outliers exist, these 
outliers barely affect the fit, slope and scatter of the rela¬ 
tion. We note that gas-depleted non-isolated galaxies are 
slightly scattered to either smaller radii or larger bary¬ 
onic masses. Since these galaxies are actually deficient 
of baryons compared to galaxies of similar stellar mass, 
we assume that the radii of these gas-depleted galaxies 
are slightly smaller than the isolated sample. 

3.3.3. Radius versus Velocity 

In Figure]^ (Bottom, Right) we show the relation be¬ 
tween radius and maximum circular velocity. This rela¬ 
tion is the least linearly correlated with a Pearson r of 
0.777, which decreases to 0.577 for the non-isolated sam¬ 
ple. The radius-velocity relation is related to the specific 
angular momentum of a galaxy. This relation is the most 
scattered of the three scaling relations. The scatter may 
be partially driven by the choice of radius, but relatively 
large scatter is predicte d by models of galaxy formation 
(Avila-Reese et al. 20081. This relation presents the most 
signihcant slope change for a scaling relation due to en¬ 
vironment, becoming shallower from 1.245 ± 0.033 for 
isolated galaxies to 1.020 ± 0.039 for non-isolated galax¬ 
ies (also see Tablej^for a list of all hts). The scatter does 
not change when we fit on the isolated sample compared 
to the non-isolated sample. The change in slope implies 
that either the effective radius decreases for non-isolated 
galaxies at fixed circular velocity, or that circular veloci¬ 
ties somehow increase. Non-isolated galaxies may be un¬ 
dergoing environmental effects that are changing the dis¬ 
tribution of luminous mass and increasing the half-light 


radius of our systems. Gas-depleted non-isolated galax¬ 
ies tend to be slightly more scattered to smaller radii 
or higher maximum circular velocities. It is not obvi¬ 
ous why the gas-depleted non-isolated galaxies would be 
scattered to higher maximum circular velocities. Since 
we also observe a similar trend with radius as a func¬ 
tion of baryonic mass, we assume that these gas-depleted 
galaxies tend to have slightly smaller effective radii than 
the isolated sample. 


3.4. Baryon Fractions 

If the atomic gas supply of galaxies is indeed depen¬ 
dent on environment, we would predict a higher mea¬ 
sured baryon fraction in our isolated low mass galaxy 
sample compared to galaxies that have had their atomic 
gas depleted by environmental processes. We define the 
baryon fraction of a galaxy as the total mass of baryons 
divided by the total mass of the galaxy (baryons plus 
dark matter). 

The total baryonic mass of a galaxy is relatively 
straightforward to estimate. However, in calculating a 
galaxy’s baryonic mass, we will always come up short of 
the actual baryonic mass due to missing contributions 
from dust, molecular gas and ionized hydrogen. The 
total galaxy mass (baryons and dark matter) is much 
more ambiguous. The total galaxy mass can either be 
estimated as the total dynamical mass directly from the 
measured galaxy dynamics (hereafter the “disk baryon 
fraction”) or from indirect estimates of the total halo 
mass (hereafter the “halo baryon fraction”). In both 
cases, we dehne baryon fraction as: 


/"baryon,disk/halo 


M, 


baryon 


M, 


dyn/halo 


( 7 ) 


In Figure we present 3 different baryon fractions: 1 
disk baryon fraction and 2 halo baryon fractions. For the 
calculations below, we assume that Vw 20 ,i is the maxi¬ 
mum circular velocity of our galaxies. 

The disk baryon fraction measures the fraction of 
baryons that have settled into the disk of t he galaxy (i.e., 
the co ndensed baryon fraction of galaxies, [Zaritsky et al.| 
2014). For disk baryon fractions, we assume that all of 
our galaxies are disk galaxies and that Vw 20 ,i has been 
measured at some radius related to the size of the HI 
disk. We calculate the total dynamical mass, which in¬ 
cludes dark matter and baryons and assumes a spherical 
distribution of mass, as 


u ^ JO N _ 17 ^ 20 , i 

-'V4dyni< Rfflj — - ^ - 


[Mq]. 


( 8 ) 


Since we do not have resolved HI rotation curves, we 
estimate the HI disk radius. We use t he empirically de¬ 
rived relation of Broeils & Rhee (1997) between Dm and 
Mhi (their Equation 3), where Dm is the HI disk diam¬ 
eter at a surface density of 1Mq/pc~^, 


log Mhi = 1.96 log -1-6.52. 


( 9 ) 


This relation is calibrated with spiral and irregu lar HI 
rich galaxie s with 10^-^ < Mhi < 10 ^°'^, (also see Wang 
et al. (2014); Hallenbeck et al. (2014)). We define the HI 
disk radius as Km = L>hi/2 . 
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(B) - Halo mass, van den Bosch 2002 
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Fig. 7.— Disk (panel A) and halo (panels B and C) baryon fractions versus circular velocity. Unbarred isolated galaxies with axis ratios 
less than 0.65 are plotted as red squares. Non-isolated galaxies with axis ratios less than 0.65 and barred isolated galaxies are plotted 
as grey dots for reference. Gas-depleted, non-isolated galaxies with axis ratios less than 0.65 are plotted as open blue dots. We plot the 
median baryon fraction as solid black horizontal lines. We plot the universal baryon fraction of 0.17 as a dotted horizontal line for reference. 
Typical uncertainties in Vw20,i /baryon,disk/halo? given purely by error propagation, are on the order of 10 km s“^ and 0.01. We also 
output the median baryon fraction in the top right of each figure. The baryon fractions in panels A and B are calculated using the circular 
velocities we measure and result in flat relations. This is most likely due to the limitations of using HI circular velocities in low mass 
galaxies where the HI gas may not probe out to Umax- The baryon fractions in panel C are measured independently of circular velocity 
but is the only panel where we And a strong relation with circular velocity. 


de Blok & Walter (2014) have shown that 50% HI ve- 
locity widths originate at roughly 5Rd, where i?D is the 
optical disk scale radius. T he near equivalence of 5Rp 
and /?Hi has been noted by de Blok & Walter (2014). 
We compare dynamical masses using these two estimates 
and find that i?Hi and 5Rd produce similar results but 
i?Hi results in a less scattered relation between circular 
velocity and baryon fraction. 

We also compare dynamical mass estimates using i?Hi 
to dynamical masses using the HI beam radius, 


Cbeam = D sin 


^beam 


2 


( 10 ) 


where D is the distance to the galaxy and Ohe&m is the 
angular size of the radio beam in 21-cm. The HI beam 
radius is the maximum radius at which we could have ob¬ 
served HI emission. Since we have not measured the ac¬ 
tual radius of the disk’s maximum circular velocity, this 
radius is the maximum radius that we could have possi¬ 
bly measured Vw 20 ,i- Again, these radii give similar re¬ 
sults, but the estimate using i?Hi has the smallest amount 
of scatter in the relation below. Therefore we chose to 
estimate only one disk baryon fraction. Note that if dark 
matter exists outside of the radius of the 21-cm obser¬ 
vations (the beam), or if the rotation velocity profile is 
increasing beyond this radius and Vw 20 ,i < Vnax, then 
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we would estimate a larger baryon fraction. 

We next calculate halo baryon fractions of our isolated 
galaxies using two diff erent estimates of viri alized halo 
mass. In their models, van den Bosch (2002|) (hereafter 
vdB02) find that the product of the square maximum cir¬ 
cular velocity and the disk scale radius provide the best 
recovery of virial mass when compared to estimates us¬ 
ing total galaxy luminosity or maximum circular velocity 
alone. We therefore calculate virial mass using scale ra¬ 
dius and maximum circular velocity via equation 15 of 
vdB02 as, 


M, 


halo,vdB02 




We assume exponential disks in our galaxies and set 
i?D = l-67reff. This is a valid approximation for most 
of our galaxies where the median Sersic index is 1.3. 

We also use the stellar mass to halo mass relation of 


Behroozi et al. (2010) (hereafter BIO). This relation be- 
tween stellar rnass and halo mass has been calibrated 
by taking advantage of the abundance matching tech¬ 
nique. This technique assumes that stellar mass is di¬ 
rectly related to halo mass with scatter less than 0.2 dex 
at all masses. For this halo mass estimate, BIO assume 
that the stellar mass function accurately maps onto halo 
mass function, even at the lowest masses where gala xies 


baryonic mass are of ten dominated by atomic gas (Pa- 


pastergis et al. 2012) and where many galaxies may be 


missed due to the surface b rightness limitations of the 
input stellar mass function (|Blanton et al. |2005b ). BIO 
abundance match SDSS stellar mass functions to halo 
mass functions with virial masses at 337 times the mean 
background density at z = 0 (the same virial mass defini¬ 
tion as vdB02). We use BIO Equation 21 with their free 
/r and k model for 0 < 2 < 1, we call this mass estimate 

Tfhalo.BlO- 

In Figure [Tj we present each of the three different esti¬ 
mates of baryon fractions. Uncertainties on Vw 20 ,i and 
/baryon are typically on the order of 10 km s“^ and 0.01 
respectively, given purely by error propagation. We list 
the median baryon fractions with rms scatter at the top 
right of each panel and as a yellow dashed horizontal line. 

The first panel (A) of Figure [^presents a disk baryon 
fraction as a function of circular velocity. The next two 
panels (B and C) present halo baryon fractions as a func¬ 
tion of circular velocity and stellar mass. For reference, 
we plot the universal baryon fraction /baryon = 0.17 in 
each panel as a horizontal dotted line. We discuss each 
panel of Figure We visually inspect the SDSS images 
of outliers in each of the three panels. 


(A) We measure a median disk baryon fraction of 0.15± 
0.17. This measurement is similar to the universal 
baryon fraction. We assume this is because we are 
measuring the dynamical mass at i?Hi rather than 
the virial mass at Rvir, so we are not measuring a 
significant mass contribution from the dark matter 
halo. Galaxies at large baryon fractions are largely 
driven by mis-measured dynamical masses due to 
bars that have gone un-identified in Galaxy Zoo. 
These bars dramatically affects our baryon fractions 
due to the way we estimate dynamical masses us¬ 


ing i?Hi and Vw 20 ,i- Galaxies at very low baryon 
fractions are mostly LSB galaxies. It also appears 
that the relation between HI mass and i?Hi begins 
to break down at small circular velocities or that low 
mass galaxies have larger baryon fractions in their 
galactic disks. Gas-depleted galaxies tend to lie at 
large baryon fractions, indicating that the dynami¬ 
cal masses can be severely underestimated when a 
galaxy is gas-depleted. 


(B) The halo baryon fr action estimate using the viri- 
“ ^ ^ (pO^ who fit 


alized halo mass of van den Bosch 
the relation between i?ciFmax Mvir. The median 
value is 0.04 ± 0.06. Outliers at large baryon frac¬ 
tions are due to barred galaxies not identified by the 
Galaxy Zoo catalog, possible mergers, and irregular 
morphologies. Galaxies at very low baryon fractions 
are mostly LSB galaxies. Galaxies at small circular 
velocities do not show the upturn in baryon frac¬ 
tion that appears in panel A, they remain constant 
at nearly all circular velocities with significant scat¬ 
ter and a slight decrease in baryon fraction at high 
masses. Gas-depleted galaxies appear unperturbed 
from the primary relation. 


(G) The halo baryon fraction estimate using t he stellar 
mass t o halo mass abundance matching of|Behroozi| 


et al. (|2010| recovers the only obvious relation be- 
tween circular velocity and baryon fraction. While 
the dynamical masses estimated in panel A and the 
halo masses estimated in panel B are calculated us¬ 
ing a power law relation with maximum circular ve¬ 
locity, the BIO stellar mass to halo mass relation 
results in somewhat of an exponential relation with 
maximum circular velocity (not shown here). This 
relation flattens out at small circular velocities (be¬ 
low ^ 70 km s“^ and asymptotes at large circu¬ 
lar velocities (above 200 km s“^ ). For galaxies 
with 80 < 4^20,1 < 100 km s“^ , the vdB02 and 
BIO baryon fractions are extremely similar because 
this is where the relations between maximum circu¬ 
lar velocity and halo mass overlap. The BIO mea¬ 
surements of baryon fractions imply that both low 
mass and high mass galaxies are severely deficient 
of baryons compared to both the cosmic value. The 
median baryon fraction using this halo mass esti¬ 
mate is 0.04 ± 0.02, which is equal to, but less scat¬ 
tered than the estimate in panel B. Gas depletion in 
galaxies does not affect this estimate of halo mass, 
so gas-depleted galaxies simply lie at smaller baryon 
fractions than the rest of the sample. 


In the cases above that estimate baryon fractions using 
halo mass, the measured baryon fraction is depleted com¬ 
pared to the universal baryon fraction. As we mentioned 
earlier, the baryon fractions of our galaxy sample must 
be somewhat depleted since we have not detected all of 
the ionized gas, molecular gas and dust in these galaxies. 
We assume that our atomic gas observations allow us to 
make reliable estimates of total baryonic mass. We dis¬ 
cuss the impact of molecular gas on our results further 
in § [4T| 

The scaling relation between maximum circ ular ve¬ 
locity and baryonic mass we measure in § 3.3 suggests 
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a weak correlation between maximum circular velocity 
and baryon fract ion: We measure Mbaryon oc 1/3.238 
r oc 1/1/1-245 jg g 3 3 jf assume that Mdynamic oc rV'^ 
then Mbaryon/-^dyn cc 1/0-435^ gbserve such a 

slope in the three panels above, i.e. this derived slope im¬ 
plies a smaller baryon fraction at low masses and a larger 
baryon fraction at high masses, but we observe a mostly 
constant baryon fraction in panels A and B above. 

Since the precise value of r in our dynamical mass 
estimate is uncertain, we could instead assume a di¬ 
rect relation between total virialized mass and maxi¬ 
mum circular velocity. We compare to the results of 
N-body simulations from Klypin et al. (2011) who find 
that Mvir oc for distinct halos (see their equation 

8 ). Given the previous paragraph, this would imply a 
roughly constant baryon fraction at all velocities, as in 
cases A and B above. 

In general, we find that baryon fraction is independent 
of environment. The lack of environmental dependence 
o n bar yon fraction is not surprising given our results in 
§ |3.3.1 We note that the scatter in baryon fraction dou- 
bles when the non-isolated galaxy sample is included. We 
find no evidence for severe baryon depletion in low mass 
galaxies when we attempt to directly measure dynamical 
or halo mass, as opposed to the baryon deficiency using 
the abundance matching result. This result motivates 
further research into the most accurate measurement of 
halo mass in the lowest galaxy mass regime. 

4. DISCUSSION AND COMPARISON TO SIMILAR STUDIES 

We discuss our results in context with similar studies, 
specifically in terms of environmental effects. We focus 
on results from the ALFALFA survey and other work 
that has inspected comparable HI observations of late- 
type galaxies at low redshift over similar galaxy mass 
ranges. 

4.1. Gas Fractions and the Gas to Stellar Mass Relation 

Our relation betw een M* and Mhi agrees well with 
Zhang et al. (2009), who use the homogeneous Hyper- 
Leda catalog a!nd the SDSS DR4 to inspect the difference 
in HI masses for 800 metal-rich and metal poor galaxies 
as a function of fixed stellar mass between IO^'^Mq and 
IO^^Mq. These authors find an offset between metal- 
poor and metal-rich galaxies with metal-poor galaxies 
showing higher HI gas fractions than their metal-rich 
sample. These results suggest that galaxies with higher 
HI gas fractions are less evolved (at least chemically), in¬ 
dicating that our isolated low mass sample at high atomic 
gas fractions are the least evolved galaxies at this mass 
scale. It may also be that these gas-rich galaxies have 
simply accreted a significant amount of low metallicity 
gas. 

The atomic gas content of SDSS DR 7 galaxies have 
be en studied ex t ensivel y usi ng the ALFALFA40 d ataset 


Huang et al. (2012a) and 


after HI2a and H12b respectively, as well as m|Maddox| 


Huang et al.| (2012b), here- 


et al. (2014). At low masses, H12a note a similar break in 
the relation of /gas to M* at M* = 10 ®Mq. They suggest 
this break is consistent with the effects of AGN feedback 
on the atomic gas content in massive galaxies. Recently, 
evidence of AGN has been found in galaxies w ith stel¬ 
lar mass as low as IO^’^^Mq (Reines et al.||20f3), so the 


effect of AGN feedback must be less pronounced in low 
mass galaxies for this to be true. 

Why such a pronounced break in the stellar mass to 
atomic gas mass relation exists at this stellar mass is not 
entirely clear and is not predicted by e xisting models of 
galaxy formation (Maddox et al.||2014). Since the break 
in the relation is independent of environment, it is most 
likely due to internal processes of galaxy formation and 
star formation. Given the fact that we detect small frac¬ 
tions of baryons in these galaxies, we could assume that 
some significant component of baryons have gone unde¬ 
tected. 

We estimate the amount of missing gas from high mass 
galaxies using the low mass end of the relation betw een 
stellar mass and atomic gas mass measured in § 3.2[ We 
find that high mass galaxies would need roughly four 
times the measured HI mass in order to account for the 
break in the power law. While it may be that high mass, 
late-type galaxies contain a significant amount of molec¬ 
ular hydrogen (e.g. the Mi lky Way has roughly a three- 
to-one Mhi to Mh 2 ratio (Draine 2011|) and M94 has a 


one-to-one Mhi to Mh^ ratio (Leroy et al. 2009)), the 
molecular gas mass in high mass spiral galaxies appears 


1 McGaugh & de Blok||1997 

Leroy et al.||2008| Saintonge 

et ai.| 
to ma 
et al. 
ized g 
Wh 
know] 

2011 ). i'he amount of molecular hydrogen re 

keup this mass deficit would be extreme (e.g., 
2008) but combining this molecular gas wii 
as may significantly reduce the discrepancy, 
ile the ionized gas content of galaxies is no 
1, the molecular gas content of galaxies ha 

quired 

Leroy 
h ion- 

t well- 
3 been 

studied extensively (e.g., IKenney & Young 

I'IN'I 

Young 

& Scoville||19911 IMcGaugh & de Blok||1997 

(Welch et al. 


M hj to Mh^ ratio of iso lated galaxies (see Fig. 15 and 16 
of Lisenfeld et al. ]MT| ). Their results imply that galax¬ 


ies inust~5e''strbngl}nnteracting in order to increase the 
ratio of molecular gas mass to neutral gas. This is most 
likely because the HI gas disk is less tightly bound than 
the molecular gas and is stripped first by environmen¬ 
tal processes. We would not expect a significant number 
of such strongly interacting galaxies in our non-isolated 
sample due to the restriction that the HI beam is not 
confused between sources. 

We estimate the molecular gas mass of our galaxies 
using Ha flux in the SDSS fiber to estimate the star for¬ 
mation rates (SFR) for our galaxies. Since the SDSS 
fiber aperture is smaller than each galaxy, we correct 
for the missing Ha flux by computing the ratio of the 
r-band flux in a fiber aperture to the total r-band flux 
of the galaxy. To estimate molecular gas mass, we as- 
su me a relation between SFR and Mh 2 (see Equation 1 


of McGaugh & Schombert (2015)). It may be that star 
formation rates are attected by environment, so this ex¬ 
ercise may reveal somethin g about environme ntal effects 
on our galaxy sample (e.g., Perea et ar]|1997). 


We find the Mh 2 to Mhi ratio of our sample lies be¬ 
tween 0.04 and 0.08 and is highly scattered but fairly 
constant at all masses. We re-run our analysis and fit¬ 
ting algorithms and find that this estimate of molecular 
gas mass does not impact our results significantly. As 
expected, we have fewer gas depleted non-isolated galax¬ 
ies and the slopes of our relations change slightly. While 
the star formation efficiency of galaxies may not be sig- 
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nificantly effected by environment (Casasola et al. ]2004t , 
the molecular gas fraction probably is. I'herefore an ern- 
pirical estimate of molecular gas based on star formation 
efficiency arguments won’t pickup on environmental ef¬ 
fects if this was true. Also, we focus our work on low 
mass galaxies and the molecular gas fractions of low mass 
galaxies are observed to be much lower given their low 
surface brightness and lower HI surface densities (Mc- 
Gaugh|[2012 Schruba et al. 2012). We plan to study the 
effect of environment on HFKs iii a future work. 

At hi gh masses our data agree well with the GASS 
survey ( Gatinella et al.|[2010 |, with their non-detections 
tracing much lower gas tractions than our sample. The 
GASS observing strategy is to observe until an atomic 
gas fraction threshold is reached. Their work can offer 
insight into non-isolated galaxy evolution at higher gas 
masses. We are focused on the low mass regime so we 
forgo a detailed comparison to the GASS survey sample, 
but we note that GASS high mass non-isolated galaxies 
are gas-depleted relative to our gas-rich sample. 


4.2. The Gas Content of Galaxies in Low Density 
Environments 

H12a has shown that galaxies with Mi, > lO^Mfr^ in 
the d enser Virgo Gluster environment (jCortese et al. 
2011 ) are gas-depleted compared to the primary, gas- 
rich A LFALFA sample (see their figure 2c). |Cortese et al. 
(2011) claim that by comparing different models of Virgo 
dwarfs, they can determine that ram-pressure stripping 
is the cause of gas depletion in Virgo spirals as opposed 
to slow starvation. What is interesting about our re¬ 
sults is that despite the ALFALFA40 sample being so 
gas-rich, we identify a small population of non-isolated 
galaxies that are noticeably gas-depleted compared to 
isolated galaxies. It would be interesting to compile a 
deep HI survey of non-isolated galaxies from the SDSS 
and probe this relation down to co mparable atomic gas 
fractions as in (jCortese et al.||2011 ). 

While some effort has been made to study the HI con¬ 
tent of void galaxies (e.g., Kreckel et al. |2012| ), our defini¬ 
tion of isolated means not all isolated galaxies are located 
in voids and not all void galaxies are isolated. Therefore 
results from studies of void galaxies may not be as use¬ 
ful as a direct comparison to our isolated galaxies. Re¬ 
gardless, Kreckel et al. (2012) find that low mass void 


galaxies with stellar mass between Mq and 10^Mq 
are not preferentially gas-rich. Our results suggest that 
depending on their definition, void galaxies may be still 
be affected by their local environment, especially at such 
low masses. 


4.3. Scaling Relations 


Our velocity-mass-size scaling relations agree well 
with similar studies of gas-rich spirals at higher masses 


(ICourteau et al. 112007 Avila-Reese et al. 

2008 

Tollerud 

et al. 2010 Hall et al. 2012). Since we iind 

that en- 


vironment can affect the atomic gas fractions of galax¬ 
ies, it would be logical that the scaling relations of these 
galaxies are also affected. While galaxy scaling relations 
become more scattered when non-isolated galaxies are in¬ 
cluded, we find that the galaxy scaling relations’ slopes 
are not significantly perturbed compared to our isolated 
gas-rich sample. 


It may not be that galaxies move perpendicular to scal- 
ing relations as the y undergo environmental processes. 
Dutton et al. (2010) suggest that observed disk galaxies 
may evolve parallel to their scaling relations, which may 
be the case for galaxy evolution due to environmental 
processes. 

For example, as a galaxy is stripped of atomic gas, 
it may not form stars at large radii where the binding 
energy is the lowest and the atomic gas is the most de¬ 
pleted. With its fuel reservoir depleted, a galaxy may 
also form fewer stars overall than a galaxy that has not 
been stripped of atomic gas. Therefore, as a galaxy de¬ 
creases in radius due to environmental processes it may 
also decrease in stellar mass, moving parallel to the scal¬ 
ing relation between effective radius and baryonic mass. 
In summary, the effect of environment on galaxies may be 
difficult to measure using the velocity-mass-size scaling 
relations. 

It may also be that we do not observe a significant 
environmental dependence in galaxy scaling relations be¬ 
cause the ALFALFA survey does not measure the scaling 
relations of gas-deficient disk galaxies due to the nature 
of a flux limited survey. In the d ensest environments , 
such as in rich group environme nts (jGatinella et al.|20I3 1 
and dense cluster environments (Kenney & Young|I989), 
the HI deficiency of galaxies increases dramatically as 
a fu nction of radius from th e X-ray center of the clus¬ 
ter (jBoselli fc Gavazzi 2006). Fernandez Lorenzo et al. 


(201^1 find evidence that massive isolated spiral galaxies 


with 10^0 < M* < IO^IMq are larger than galaxies in 
denser environments at fixed mass. 

When SO galaxies are compared to the Tully-Fisher re¬ 
lation of spirals, an offset and an increase in scatter are 


observed (Neistein et al.|I999 

Hinz et al. 12003 

Bedregall 

et al. 

2006 Davis et ai.||20II 

Williams et al. 

2010 |. it 


turbed, compact or merging they do not follow the Tully- 
Fisher relation and the circular velocity must be added 


to the non-circular motions of the HI gas (Weiner et al. 
2006 Kassin et al.||2007 2014). Taking into account the 


velocity dispersion and concentration of a galaxy, disk 
and spheroidal galaxi es follow similar scaling relations 
(Gatinella et al. 2012). If spheroidal galaxies are trans- 
lormed by their environment, then galaxies in groups and 
clusters should systematically vary fr om the scaling rela- 
tions presented in this work. Indeed, Gortes et al. (2008) 
claim that environmental effects have decreased HI line 
widths in Virgo cluster galaxies. However, the most no¬ 
ticeable environmental effect on the galaxy scaling rela¬ 
tions is that gas-depleted non-isolated galaxies tend to be 
scattered to smaller effective radii than isolated galaxies 
- indicating again that we are not sensitive to environ¬ 
mental effects on circular velocity or that Vw 20 ,i is also 
encapsulating noncircular motion. 

4.4. Baryon Fractions 

We find that baryon fraction is independent of circu¬ 
lar velocity in all cases except Mhaio.Bio- The baryon 
fraction calculated using Mhaio.Bio shows a significant 
dependence on circular velocity because this estimate is 
based on stellar-t o-halo mass abundance m atching. This 
result is similar to Papastergis et al. (2012), who perform 
abundance matching with total baryonic masses using 
ALFALFA40 HI masses. These authors find that the 
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depletion of baryons at low halo masses is only mildly 
alleviated when total baryon mass is used in the abun¬ 
dance matching technique. It has been suggested that 
LSB galaxies may significantly change abundance match¬ 
ing results by changing the stellar mass function that 


goes into matching the halo mass function (Ferrero et al. 


20121 . 

Dur other baryon fraction estimates are at odds with 
the abundance matching technique which suggest baryon 
fractions are dependent on circular velocity and are much 
more significantly depleted of baryons. The flat relations 
we measure between circular velocity and baryon frac¬ 
tions at low circular velocities suggests that our HI veloc¬ 
ity width measurement is not probi ng deep into the dar k 
matter halo of low mass galaxies (Ferrero et al. 2012). 
This could be confirmed by resolved HI synthesis maps. 
If the HI rotation curves do not reach the maximum cir¬ 
cular velocity of the dar k matter halo (e.g., s everal exam¬ 
ples are available from |Swaters et al.||2009|), our Vw 20 ,i 
measurement would underestimate tiaio mass. This ef- 
fe ct would lead to a flat tening of the relation. 


Blanton et al. (2008) find a weak relation between 
baryon traction as a fu nction of luminosity at low masses. 


Zaritsky et al. (2014) estimate the disk baryon fraction 
tor a simi lar sample of gas-rich galaxies using the mass es- 
|2 001 1 and find a baryon fraction 


timate of Bullock et al. 


of 0.07. McGaugh et al. 1 2009) calculate the fraction of 
detected baryons to range from 0.83 for clusters to 0.0003 
for the lowest mass dwarf galaxies. Our results do not in¬ 
dicate such a significant depletion of baryons in isolated 
galaxies with stellar mass 10^ < < 10® M©, nor do 

we note a significant trend towards low baryon fractions 
except in the BIO abundance matching estimate. 


5. SUMMARY 

We present new measurements of HI gas masses for 
148 isolated galaxies with stellar masses between 10^ and 
IO®-®M0. We apply the same isolation criterion to our 
low mass sample as G12 and extend this isolation cri¬ 
teria to galaxies at higher masses. We re-measure the 
HI masses and velocity widths of galaxies observed with 
the ALFALFA survey in order to create a large, homoge¬ 
neous data set. Our HI data set represents the largest of 
its kind with such a rigorously defined isolation criteria. 

Using these data, we study the baryon content of galax¬ 
ies in isolated environments, specifically at low masses. 
We inspect the relation between galaxy stellar mass and 
atomic gas mass in isolated environments versus the rela¬ 
tion in non-isolated galaxies (i.e. denser environments). 
We also examine the relations between velocity, mass and 
size. Finally, we determine the disk baryon fractions and 
two estimates of halo baryon fractions for our galaxies. 
Our primary results are as follows: 


• We either place all isolated low mass galaxies with 
stellar masses between 10^ and I0 ® ®Mq at high 
atomic gas fractions or we cannot rule out high 
atomic gas fractions via upper limits. Isolated 
galaxies have a median atomic gas fraction /gas = 
0.81 ±0.13 and atomic gas fractions always greater 
than 0.3. 


• We use a largely flux-limited sample of non-isolated 
low mass galaxies to derive a similar median atomic 


gas fraction as isolated galaxies. We measure non¬ 
isolated atomic gas fractions as low as 0.15, which 
is limited by the Hl-rich galaxies of the ALFALFA 
survey. 

• We observe a break in the relation between stellar 

mass and atomic gas mass at 10® Both iso¬ 

lated and non-isolated galaxies follow this relation. 

• Only non-isolated galaxies are depleted in atomic 
gas by more than —1 dex compared to the atomic 
gas mass predicted by the stellar mass to atomic 
gas mass relation. 

• Environment does not significantly affect the linear 
correlations between velocity, mass and size. Gas- 
depleted non-isolated galaxies tend to be scattered 
to smaller effective radii. 

• Isolated galaxies follow a steeper relation between 
velocity and radius than non-isolated galaxies. 

• The baryon fraction of isolated galaxies is inde¬ 
pendent of environment: we measure a median 
baryon to dynamical mass fraction of /baryon,disk = 
0.15 ± 0.17 and median baryon to halo mass frac¬ 
tions of /baryon,halo = 0.04±0.06 using the results of 
semi-analytic models and /baryon,halo = 0.04 ±0.02 
using abundance matching. 

The above results improve our understanding of the 
baryon content of low mass galaxies and offer us some 
insight into the effect of environment in this stellar mass 
regime. Extending our analysis to galaxies below our 
stellar mass limit of Mq requires deeper spectroscopic 
surveys (e.g., DESI) and deeper all-sky HI surveys (e.g., 
SKA). The flat relation we measure between circular ve¬ 
locity and baryon fraction motivates resolved, follow-up 
HI synthesis maps of our low mass galaxy sample in or¬ 
der to determine the true maximum circular velocity of 
our systems and the radius where we measure this ve¬ 
locity. We plan to study the spatial distribution and the 
kinematics of baryons in our sample using integral field 
unit observations. 
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TABLE 1 

Sample of Our High Quality HI Measurements and Environments for Isolated Low Mass Galaxies 


NSA ID 

(1) 

a (J2000.0) 
[h: m: s] 
(2) 

5 (J2000.0) 

(3) 

Distance 

[Mpc] 

(4) 

Mr 

[mag] 

(5) 

VEso 

[km s ^ ] 
(6) 

S 21 

[Jy km s“^[ 

(7) 

Mhi 

[1O®M0[ 

(8) 

M, 

[1O®M0[ 

(9) 

^host 

[Mpc] 

(10) 

162 

09:58:30.23 

+00:02:42.72 

30.23+1.89 

-13.34 

50+5 

0.49+0.03 

0.105+0.029 

0.016+0.007 

2.15 

1443 

12:23:09.98 

+00:25:37.76 

35.36+2.00 

-13.61 

32+28 

0.14+0.03 

0.041+0.019 

0.043+0.020 

1.69 

5109 

22:30:36.80 

-00:06:37.01 

21.26+1.86 

-13.40 

55+3 

0.63+0.02 

0.067+0.024 

0.008+0.004 

6.28 

5378 

23:06:39.16 

+00:34:41.02 

68.13+2.39 

-15.86 

104+7 

0.95+0.03 

1.039+0.166 

0.273+0.126 

1.77 

5424 

23:14:19.43 

+00:10:59.36 

60.30+1.75 

-15.67 

108+7 

0.47+0.02 

0.405+0.058 

0.612+0.282 

2.43 

5449 

23:28:12.31 

-01:03:46.37 

35.30+2.14 

-13.90 

64+2 

0.78+0.02 

0.229+0.059 

0.020+0.009 

1.51 

5720 

00:13:44.02 

+00:22:18.46 

56.50+1.75 

-14.63 

29+2 

0.45+0.02 

0.337+0.049 

0.097+0.045 

1.87 

6186 

00:57:56.60 

+00:52:08.56 

30.17+2.86 

-13.91 

58+2 

0.70+0.02 

0.150+0.059 

0.051+0.024 

3.26 

6368 

01:09:07.85 

+01:07:15.89 

15.97+1.64 

-13.03 

40+1 

1.02+0.02 

0.061+0.026 

0.015+0.007 

1.69 

6412 

01:13:45.98 

+00:56:44.13 

104.68+2.14 

-16.92 

131+21 

0.30+0.02 

0.766+0.135 

0.788+0.363 

3.63 

6748 

01:26:05.01 

+00:19:00.68 

24.85+2.00 

-14.26 

117+1 

3.17+0.03 

0.461+0.152 

0.088+0.040 

1.63 

6949 

01:55:59.76 

-00:11:08.02 

51.69+1.75 

-15.74 

90+16 

0.33+0.02 

0.208+0.039 

0.323+0.149 

3.48 

7167 

02:18:08.02 

+00:45:29.39 

39.84+2.32 

-14.85 

116+1 

0.99+0.02 

0.371+0.091 

0.139+0.064 

3.57 

8331 

00:49:22.10 

+15:31:33.00 

76.64+2.89 

-15.81 

126+4 

0.77+0.02 

1.059+0.178 

0.400+0.184 

2.18 

8333 

00:56:42.19 

+14:12:45.12 

156.79+2.14 

-17.68 

88+37 

0.67+0.10 

3.883+1.242 

2.932+1.350 

3.21 

Note. — 

Properites of 15 i 

solated, low mass NSA galaxies that 

have been 

observed in HI. 

Columns are (1) 

the NASA Sloan 

Atlas vO.1.2 (NSA) 

identifier 


for all galaxies which is unique to this iteration of the NSA catalog, (2) right ascension, (3) declination, (4) heliocentric distance, (5) absolute r-magnitude, 
(6) width of the 21-cm line at 50 per cent flux, (7) integrated flux of the 21*cm line, (8) HI mass, (9) stellar mass and (10) distance to the nearest massive 
host. This table is published in its entirety in the electronic edition of the Astronomical Journal, A portion is shown here for guidance regarding its form and 
content. 


TABLE 2 

HI Observation Details 


Run ID 

Avg. Integration Time [s] 

N Reduced 

N High Quality 

N Non-Detections 

N Bad Baselines 

GBT2008 

802 

184 

141 

43 

12 

A02005 

796 

32 

32 

0 

3 

AO2013 

459 

83 

80 

3 

8 

AO2014a 

499 

51 

47 

4 

0 

AO2014b 

571 

21 

17 

4 

0 


Note. — Observing run identifier, average integration time, number of reduced spectra, number of high quality spectra, number of 
non-detections, number of bad baselines for all of our HI observations. Note that GBT has a smaller dish diameter than Arecibo, therefore 
integration times are not directly comparable. 


TABLE 3 

Scaling Relations with log y = a log x + /3 ± e 


sample 

X 

y 

a (Slope) 

(3 (Constant) 

e (Scatter in logy) 

Pearson r 

Figure 

isolated only 

M* < 1O®-®M0 

Afgas 

1.052+ 0.058 

0.236+ 0.476 

0.285+ 0.019 

N/A 


5 


isolated only 

M* > 1O«-®M0 

Afgas 

0.461+ 0.011 

5.329+ 0.112 

0.221+ 0.006 

N/A 


S 


isolated only 


TefI 

0.406+ 0.007 

-3.456+ 0.074 

0.159+ 0.004 

0.880 


6 


non-isolated only 

-^^baryon 

Teff 

0.390+ 0.009 

-3.293+ 0.094 

0.179+ 0.004 

0.764 

N 

7 

A 

all 


Teff 

0.400+ 0.006 

-3.386+ 0.059 

0.171+ 0.003 

0.826 

N/ 

A 

isolated only 

Vw20,i 

TefI 

1.245+ 0.033 

-2.005+ 0.069 

0.209+ 0.005 

0.777 


6 

n 


non-isolated only 

'^W20,i 

TefI 

1.020+ 0.039 

-1.522+ 0.083 

0.225+ 0.005 

0.577 

N 


A 

all 

Vw20,i 


1.147+ 0.026 

-1.793+ 0.054 

0.219+ 0.004 

0.683 

N/ 

A 

isolated only 

-^^baryon 

Ew20,i 

0.277+ 0.004 

-0.672+ 0.041 

0.075+ 0.002 

0.916 


6 

n 


non-isolated only 

-^^baryon 

Vw 20 ,i 

0.266+ 0.005 

-0.565+ 0.050 

0.080+ 0.002 

0.831 

N 


A 

all 

^^baryon 

yW20,i 

0.272+ 0.003 

-0.624+ 0.031 

0.078+ 0.002 

0.878 

N/A 

isolated only 

bv20,i 

-^^baryon 

3.238+ 0.046 

3.210+ 0.097 

0.250+ 0.008 

0.916 

N/A 

non-isolated only 

Vw 20 ,i 

-^^baryon 

2.949+ 0.053 

3.834+ 0.112 

0.261+ 0.007 

0.831 

N/A 

all 

'^W20,i 

-^^baryon 

3.113+ 0.035 

3.482+ 0.074 

0.258+ 0.006 

0.878 

N/A 

Note. — All linear r 

lations fit in this pape 

. Sample, do 

main, range, slope. 

constant, scatter. 

Pearson r and Figure reference for each fit 

For 

all 

mass-velocity-radius scaling relations, we have 

excluded galaxies with bars and 

with axis ratios gre 

ater than 0.65 from the galaxy scaling rela1 

ion 

fits. 
















